Ideas:

1 Transects

Plant Flowers Date lon lat ele Month Year julian
Glossoloma oblongicalyx 4 2015-10-19 -78.59093 0.130838 2270 October 2015 292
Gasteranthus quitensis 2 2016-10-17 -78.59770 0.120070 1940 October 2016 291
Kohleria affinis 1 2016-12-13 -78.59534 0.126746 2110 December 2016 348
Columnea ciliata 3 2014-02-27 -78.59934 0.116682 1960 February 2014 58
Columnea medicinalis 1 2014-04-23 -78.59372 0.128700 2130 April 2014 113
Drymonia teuscheri 3 2016-07-28 -78.59245 0.129393 2200 July 2016 210

2 Interactions

3 Phylogeny

4 Traits

4.0.1 Total Flowers

4.1 Peak date

As range

4.2 Species elevation ranges

4.3 Flowering Data Matrix

5 Baseline model

Equal probability of flowering at anytime.

## sink("model/threshold_baseline.jags")
## cat("
##     model {
##     
##     for (x in 1:Nobs){
## 
##     #Observation of a flowering plant
##     Y[x] ~ dbern(p[x])
##     logit(p[x]) <- alpha[Plant[x]] 
##     
##     #Residuals
##     discrepancy[x] <- abs(Y[x] - p[x])
##     
##     #Assess Model Fit
##     Ynew[x] ~ dbern(p[x])
##     discrepancy.new[x]<-abs(Ynew[x] - p[x])
##     }
##     
##     
##     #Sum discrepancy
##     fit<-sum(discrepancy)/Nobs
##     fitnew<-sum(discrepancy.new)/Nobs
##     
##     #Prediction
##     
##     for(x in 1:Npreds){
##     #predict value
##     
##     #Observation - probability of flowering
##     prediction[x] ~ dbern(p_new[x])
##     logit(p_new[x])<-alpha[PredPlant[x]]
##     
##     #predictive error
##     pred_error[x] <- abs(Ypred[x] - prediction[x])
##     }
##     
##     
##     #Predictive Error
##     fitpred<-sum(pred_error)/Npreds
##     
##     #Priors
##     
##     #Species level priors
##     
##     for (j in 1:Plants){
##     
##     #Intercept flowering probability
##     alpha[j] ~ dnorm(0,0.386)
##     
##     } 
## 
##     }
##     ",fill=TRUE)
## 
## sink()
## Compiling model graph
##    Resolving undeclared variables
##    Allocating nodes
## Graph information:
##    Observed stochastic nodes: 2324
##    Unobserved stochastic nodes: 2926
##    Total graph size: 14658
## 
## Initializing model

6 Get Chains

6.0.1 Evaluate convergence

6.0.2 Posterior estimates

7 Phylogeny

7.1 Attraction

## sink("model/threshold_attraction.jags")
## cat("
##     model {
##     
##     for (x in 1:Nobs){
##     
##     #Observation of a flowering plant
##     Y[x] ~ dbern(p[x])
##     logit(p[x]) <- alpha[Plant[x]] + e[Plant[x],Month[x]]
##     
##     #Residuals
##     discrepancy[x] <- abs(Y[x] - p[x])
##     
##     #Assess Model Fit
##     Ynew[x] ~ dbern(p[x])
##     discrepancy.new[x]<-abs(Ynew[x] - p[x])
##     }
##     
##     #Sum discrepancy
##     fit<-sum(discrepancy)/Nobs
##     fitnew<-sum(discrepancy.new)/Nobs
##     
##     #Prediction
##     
##     for(x in 1:Npreds){
##     #predict value
##     
##     #Observation - probability of flowering
##     prediction[x] ~ dbern(p_new[x])
##     logit(p_new[x])<- alpha[NewPlant[x]] + e[NewPlant[x],NewMonth[x]]
##     
##     #predictive error
##     pred_error[x] <- abs(Ypred[x] - prediction[x])
##     }
##     
##     #Predictive Error
##     fitpred<-sum(pred_error)/Npreds
##     
##     #########################
##     #autocorrelation in error
##     #########################
##     
##     #For each of observation
##     for(y in 1:Months){
##     e[1:Plants,y] ~ dmnorm(zeros,tauC[,])
##     }
##     
##     ##covariance among similiar species
##     for(i in 1:Plants){
##     for(j in 1:Plants){
##     C[i,j] = exp(-lambda_cov * D[i,j])
##     }
##     }
##     
##     ## Covert variance to precision for each parameter, allow omega to shrink to identity matrix
##     vCov = omega*C[,] + (1-omega) * I
##     tauC=inverse(vCov*gamma)
##     
##     #Priors
##     
##     #Species level priors
##     
##     for (j in 1:Plants){
##     
##     #Intercept flowering probability
##     alpha[j] ~ dnorm(0,0.386)
##     
##     } 
##     #Autocorrelation priors
##     gamma  = 100
##     
##     #Strength of covariance decay
##     lambda_cov = 5
##     omega = 1
##     }
##     ",fill=TRUE)
## 
## sink()
## Compiling model graph
##    Resolving undeclared variables
##    Allocating nodes
## Graph information:
##    Observed stochastic nodes: 2324
##    Unobserved stochastic nodes: 2938
##    Total graph size: 19058
## 
## Initializing model

8 Get Chains

8.0.1 Evaluate convergence

8.0.2 Posterior estimates

8.1 Effect of phylogenetic attraction

8.2 Repulsion

## sink("model/threshold_repulsion.jags")
## cat("
##     model {
##     
##     for (x in 1:Nobs){
##     
##     #Observation of a flowering plant
##     Y[x] ~ dbern(p[x])
##     logit(p[x]) <-  alpha[Plant[x]] + e[Plant[x],Month[x]]
##     
##     #Residuals
##     discrepancy[x] <- abs(Y[x] - p[x])
##     
##     #Assess Model Fit
##     Ynew[x] ~ dbern(p[x])
##     discrepancy.new[x]<-abs(Ynew[x] - p[x])
##     }
##     
##     
##     #Sum discrepancy
##     fit<-sum(discrepancy)/Nobs
##     fitnew<-sum(discrepancy.new)/Nobs
##     
##     #Prediction
##     
##     for(x in 1:Npreds){
##     #predict value
##     
##     #Observation - probability of flowering
##     prediction[x] ~ dbern(p_new[x])
##     logit(p_new[x])<- alpha[NewPlant[x]] +  e[NewPlant[x],NewMonth[x]]
##     
##     #predictive error
##     pred_error[x] <- abs(Ypred[x] - prediction[x])
##     }
##     
##     #Predictive Error
##     fitpred<-sum(pred_error)/Npreds
##     
##     #########################
##     #autocorrelation in error
##     #########################
##     
##     #For each of observation
##     for(y in 1:Months){
##     e[1:Plants,y] ~ dmnorm(zeros,tauC[,])
##     }
##     
##     ##covariance among similiar species
##     for(i in 1:Plants){
##     for(j in 1:Plants){
##     C[i,j] = exp(-lambda_cov * D[i,j])
##     }
##     }
##     
##     ## Covert variance to precision for each parameter, allow omega to shrink to identity matrix
##     vCov = omega*C[,] + (1-omega) * I
##     tauC=vCov*gamma
##     
##     #Priors
##     
##     #Species level priors
##     
##     for (j in 1:Plants){
##     
##     #Intercept flowering probability
##     alpha[j] ~ dnorm(0,0.386)
##     
##     } 
##     #Autocorrelation priors
##     gamma = 0.1
##     
##     #Strength of covariance decay
##     lambda_cov = 1
##     omega = 1
##     }
##     ",fill=TRUE)
## 
## sink()
## Compiling model graph
##    Resolving undeclared variables
##    Allocating nodes
## Graph information:
##    Observed stochastic nodes: 2324
##    Unobserved stochastic nodes: 2938
##    Total graph size: 19056
## 
## Initializing model

9 Get Chains

9.0.1 Evaluate convergence

9.0.2 Posterior estimates

9.1 Effect of phylogenetic repulsion

10 Traits

10.1 Trait Attraction

## sink("model/threshold_attraction.jags")
## cat("
##     model {
##     
##     for (x in 1:Nobs){
##     
##     #Observation of a flowering plant
##     Y[x] ~ dbern(p[x])
##     logit(p[x]) <- alpha[Plant[x]] + e[Plant[x],Month[x]]
##     
##     #Residuals
##     discrepancy[x] <- abs(Y[x] - p[x])
##     
##     #Assess Model Fit
##     Ynew[x] ~ dbern(p[x])
##     discrepancy.new[x]<-abs(Ynew[x] - p[x])
##     }
##     
##     #Sum discrepancy
##     fit<-sum(discrepancy)/Nobs
##     fitnew<-sum(discrepancy.new)/Nobs
##     
##     #Prediction
##     
##     for(x in 1:Npreds){
##     #predict value
##     
##     #Observation - probability of flowering
##     prediction[x] ~ dbern(p_new[x])
##     logit(p_new[x])<- alpha[NewPlant[x]] + e[NewPlant[x],NewMonth[x]]
##     
##     #predictive error
##     pred_error[x] <- abs(Ypred[x] - prediction[x])
##     }
##     
##     #Predictive Error
##     fitpred<-sum(pred_error)/Npreds
##     
##     #########################
##     #autocorrelation in error
##     #########################
##     
##     #For each of observation
##     for(y in 1:Months){
##     e[1:Plants,y] ~ dmnorm(zeros,tauC[,])
##     }
##     
##     ##covariance among similiar species
##     for(i in 1:Plants){
##     for(j in 1:Plants){
##     C[i,j] = exp(-lambda_cov * D[i,j])
##     }
##     }
##     
##     ## Covert variance to precision for each parameter, allow omega to shrink to identity matrix
##     vCov = omega*C[,] + (1-omega) * I
##     tauC=inverse(vCov*gamma)
##     
##     #Priors
##     
##     #Species level priors
##     
##     for (j in 1:Plants){
##     
##     #Intercept flowering probability
##     alpha[j] ~ dnorm(0,0.386)
##     
##     } 
##     #Autocorrelation priors
##     gamma  = 100
##     
##     #Strength of covariance decay
##     lambda_cov = 5
##     omega = 1
##     }
##     ",fill=TRUE)
## 
## sink()
## Compiling model graph
##    Resolving undeclared variables
##    Allocating nodes
## Graph information:
##    Observed stochastic nodes: 2324
##    Unobserved stochastic nodes: 2938
##    Total graph size: 19214
## 
## Initializing model

11 Get Chains

11.0.1 Evaluate convergence

11.0.2 Posterior estimates

11.1 Effect of trait attraction

11.2 Repulsion

## sink("model/threshold_repulsion.jags")
## cat("
##     model {
##     
##     for (x in 1:Nobs){
##     
##     #Observation of a flowering plant
##     Y[x] ~ dbern(p[x])
##     logit(p[x]) <-  alpha[Plant[x]] + e[Plant[x],Month[x]]
##     
##     #Residuals
##     discrepancy[x] <- abs(Y[x] - p[x])
##     
##     #Assess Model Fit
##     Ynew[x] ~ dbern(p[x])
##     discrepancy.new[x]<-abs(Ynew[x] - p[x])
##     }
##     
##     
##     #Sum discrepancy
##     fit<-sum(discrepancy)/Nobs
##     fitnew<-sum(discrepancy.new)/Nobs
##     
##     #Prediction
##     
##     for(x in 1:Npreds){
##     #predict value
##     
##     #Observation - probability of flowering
##     prediction[x] ~ dbern(p_new[x])
##     logit(p_new[x])<- alpha[NewPlant[x]] +  e[NewPlant[x],NewMonth[x]]
##     
##     #predictive error
##     pred_error[x] <- abs(Ypred[x] - prediction[x])
##     }
##     
##     #Predictive Error
##     fitpred<-sum(pred_error)/Npreds
##     
##     #########################
##     #autocorrelation in error
##     #########################
##     
##     #For each of observation
##     for(y in 1:Months){
##     e[1:Plants,y] ~ dmnorm(zeros,tauC[,])
##     }
##     
##     ##covariance among similiar species
##     for(i in 1:Plants){
##     for(j in 1:Plants){
##     C[i,j] = exp(-lambda_cov * D[i,j])
##     }
##     }
##     
##     ## Covert variance to precision for each parameter, allow omega to shrink to identity matrix
##     vCov = omega*C[,] + (1-omega) * I
##     tauC=vCov*gamma
##     
##     #Priors
##     
##     #Species level priors
##     
##     for (j in 1:Plants){
##     
##     #Intercept flowering probability
##     alpha[j] ~ dnorm(0,0.386)
##     
##     } 
##     #Autocorrelation priors
##     gamma = 0.1
##     
##     #Strength of covariance decay
##     lambda_cov = 1
##     omega = 1
##     }
##     ",fill=TRUE)
## 
## sink()
## Compiling model graph
##    Resolving undeclared variables
##    Allocating nodes
## Graph information:
##    Observed stochastic nodes: 2324
##    Unobserved stochastic nodes: 2938
##    Total graph size: 19212
## 
## Initializing model

12 Get Chains

12.0.1 Evaluate convergence

12.0.2 Posterior estimates

12.1 Effect of trait repulsion

13 Interaction

13.1 Attraction

## sink("model/threshold_attraction.jags")
## cat("
##     model {
##     
##     for (x in 1:Nobs){
##     
##     #Observation of a flowering plant
##     Y[x] ~ dbern(p[x])
##     logit(p[x]) <- alpha[Plant[x]] + e[Plant[x],Month[x]]
##     
##     #Residuals
##     discrepancy[x] <- abs(Y[x] - p[x])
##     
##     #Assess Model Fit
##     Ynew[x] ~ dbern(p[x])
##     discrepancy.new[x]<-abs(Ynew[x] - p[x])
##     }
##     
##     #Sum discrepancy
##     fit<-sum(discrepancy)/Nobs
##     fitnew<-sum(discrepancy.new)/Nobs
##     
##     #Prediction
##     
##     for(x in 1:Npreds){
##     #predict value
##     
##     #Observation - probability of flowering
##     prediction[x] ~ dbern(p_new[x])
##     logit(p_new[x])<- alpha[NewPlant[x]] + e[NewPlant[x],NewMonth[x]]
##     
##     #predictive error
##     pred_error[x] <- abs(Ypred[x] - prediction[x])
##     }
##     
##     #Predictive Error
##     fitpred<-sum(pred_error)/Npreds
##     
##     #########################
##     #autocorrelation in error
##     #########################
##     
##     #For each of observation
##     for(y in 1:Months){
##     e[1:Plants,y] ~ dmnorm(zeros,tauC[,])
##     }
##     
##     ##covariance among similiar species
##     for(i in 1:Plants){
##     for(j in 1:Plants){
##     C[i,j] = exp(-lambda_cov * D[i,j])
##     }
##     }
##     
##     ## Covert variance to precision for each parameter, allow omega to shrink to identity matrix
##     vCov = omega*C[,] + (1-omega) * I
##     tauC=inverse(vCov*gamma)
##     
##     #Priors
##     
##     #Species level priors
##     
##     for (j in 1:Plants){
##     
##     #Intercept flowering probability
##     alpha[j] ~ dnorm(0,0.386)
##     
##     } 
##     #Autocorrelation priors
##     gamma  = 100
##     
##     #Strength of covariance decay
##     lambda_cov = 5
##     omega = 1
##     }
##     ",fill=TRUE)
## 
## sink()
## Compiling model graph
##    Resolving undeclared variables
##    Allocating nodes
## Graph information:
##    Observed stochastic nodes: 2324
##    Unobserved stochastic nodes: 2938
##    Total graph size: 19214
## 
## Initializing model

13.1.1 Evaluate convergence

13.1.2 Posterior estimates

13.2 Effect of interaction attraction

13.3 Repulsion

## sink("model/threshold_repulsion.jags")
## cat("
##     model {
##     
##     for (x in 1:Nobs){
##     
##     #Observation of a flowering plant
##     Y[x] ~ dbern(p[x])
##     logit(p[x]) <-  alpha[Plant[x]] + e[Plant[x],Month[x]]
##     
##     #Residuals
##     discrepancy[x] <- abs(Y[x] - p[x])
##     
##     #Assess Model Fit
##     Ynew[x] ~ dbern(p[x])
##     discrepancy.new[x]<-abs(Ynew[x] - p[x])
##     }
##     
##     
##     #Sum discrepancy
##     fit<-sum(discrepancy)/Nobs
##     fitnew<-sum(discrepancy.new)/Nobs
##     
##     #Prediction
##     
##     for(x in 1:Npreds){
##     #predict value
##     
##     #Observation - probability of flowering
##     prediction[x] ~ dbern(p_new[x])
##     logit(p_new[x])<- alpha[NewPlant[x]] +  e[NewPlant[x],NewMonth[x]]
##     
##     #predictive error
##     pred_error[x] <- abs(Ypred[x] - prediction[x])
##     }
##     
##     #Predictive Error
##     fitpred<-sum(pred_error)/Npreds
##     
##     #########################
##     #autocorrelation in error
##     #########################
##     
##     #For each of observation
##     for(y in 1:Months){
##     e[1:Plants,y] ~ dmnorm(zeros,tauC[,])
##     }
##     
##     ##covariance among similiar species
##     for(i in 1:Plants){
##     for(j in 1:Plants){
##     C[i,j] = exp(-lambda_cov * D[i,j])
##     }
##     }
##     
##     ## Covert variance to precision for each parameter, allow omega to shrink to identity matrix
##     vCov = omega*C[,] + (1-omega) * I
##     tauC=vCov*gamma
##     
##     #Priors
##     
##     #Species level priors
##     
##     for (j in 1:Plants){
##     
##     #Intercept flowering probability
##     alpha[j] ~ dnorm(0,0.386)
##     
##     } 
##     #Autocorrelation priors
##     gamma = 0.1
##     
##     #Strength of covariance decay
##     lambda_cov = 1
##     omega = 1
##     }
##     ",fill=TRUE)
## 
## sink()
## Compiling model graph
##    Resolving undeclared variables
##    Allocating nodes
## Graph information:
##    Observed stochastic nodes: 2324
##    Unobserved stochastic nodes: 2938
##    Total graph size: 19212
## 
## Initializing model

14 Get Chains

14.0.1 Evaluate convergence

14.0.2 Posterior estimates

14.1 Effect of interaction repulsion

15 Model Comparison

15.0.1 E: The effect of autocorrelation

16 One example

Glossoloma purpureum

17 By site

17.1 1300m - 1500m

17.2 1500-1700m

17.3 1700-1900

17.4 1900-2100m

17.5 2100-2300m

17.6 2300m - 2500m

17.7 Alpha: Species specific flowering rate

17.8 Omega: The magnitude of the effect of autocorrelation on mean flowering occurrence

17.9 Gamma: The variance of the effect of autocorrelation on mean flowering occurrence

17.10 Lambda: The decay in autocorrelation effect

17.11 Decay in autocorrelation effect

18 Model Fit

18.1 Bayesian pvalue

## # A tibble: 7 x 2
##   Model                        p
##   <chr>                    <dbl>
## 1 baseline                0.65  
## 2 interaction_attraction  0.0275
## 3 interaction_repulsion   0.0775
## 4 phylogenetic_attraction 0.025 
## 5 phylogenetic_repulsion  0.04  
## 6 trait_attraction        0.035 
## 7 trait_repulsion         0.0525

Without baseline 

## # A tibble: 6 x 2
##   Model                        p
##   <chr>                    <dbl>
## 1 interaction_attraction  0.0275
## 2 interaction_repulsion   0.0775
## 3 phylogenetic_attraction 0.025 
## 4 phylogenetic_repulsion  0.04  
## 5 trait_attraction        0.035 
## 6 trait_repulsion         0.0525

18.2 Model Fit

Model mean lower upper
interaction_repulsion 0.2244776 0.2168625 0.2313833
trait_repulsion 0.2239220 0.2167097 0.2310252
phylogenetic_repulsion 0.2235663 0.2165512 0.2309017
interaction_attraction 0.2204553 0.2132964 0.2272275
phylogenetic_attraction 0.2202904 0.2132722 0.2283090
trait_attraction 0.2202308 0.2132088 0.2291204

18.2.1 Without baseline

18.3 By Species

18.3.1 Without baseline

18.3.2 Zoom in

19 Prediction

19.0.1 Tables

Model mean lower upper
baseline 0.2485204 0.2244898 0.2738946
interaction_repulsion 0.2298469 0.2074830 0.2551020
trait_repulsion 0.2293070 0.2057823 0.2568027
phylogenetic_repulsion 0.2281420 0.2023810 0.2517007
phylogenetic_attraction 0.2248384 0.2023810 0.2517007
interaction_attraction 0.2233036 0.2006803 0.2448980
trait_attraction 0.2228401 0.1972789 0.2500000

19.0.2 Flowering Rates among models

Plant Month Observed baseline interaction_attraction interaction_repulsion phylogenetic_attraction phylogenetic_repulsion trait_attraction trait_repulsion
Besleria solanoides Apr 13.3 17.4 13.0 13.4 14.8 13.7 12.6 14.1
Besleria solanoides Aug 8.3 16.6 7.8 10.0 9.3 9.6 7.9 9.2
Besleria solanoides Dec 12.5 15.3 12.0 18.2 12.4 16.4 14.0 13.6
Besleria solanoides Feb 0.0 16.8 1.8 1.5 1.3 0.9 1.2 1.1
Besleria solanoides Jan 0.0 17.3 0.9 1.9 1.0 1.7 1.7 1.2
Besleria solanoides Jul 6.7 16.7 8.2 9.0 5.9 7.7 6.4 9.3
Besleria solanoides Jun 10.0 16.8 10.1 10.1 9.4 10.5 9.4 10.9
Besleria solanoides Mar 0.0 16.2 1.4 3.0 1.2 1.1 0.8 2.1
Besleria solanoides May 11.1 16.7 10.6 10.7 10.6 12.0 10.9 12.2
Besleria solanoides Nov 33.3 15.9 32.1 36.2 33.6 35.0 33.0 36.2
Besleria solanoides Oct 57.1 16.6 55.8 57.0 57.2 58.3 56.2 57.3
Besleria solanoides Sep 46.7 16.8 45.9 47.1 46.7 45.7 45.9 46.9
Columnea ciliata Apr 46.7 17.2 46.8 44.8 47.4 44.9 48.4 46.5
Columnea ciliata Aug 0.0 17.6 0.9 2.4 1.5 2.7 1.0 0.9
Columnea ciliata Dec 0.0 18.0 0.9 5.9 1.3 2.9 0.5 6.3
Columnea ciliata Feb 43.8 18.5 43.3 43.8 42.6 42.7 44.9 43.2
Columnea ciliata Jan 25.0 18.2 23.8 24.2 26.0 26.3 24.4 24.3
Columnea ciliata Jul 6.7 18.4 6.2 7.8 6.1 7.5 6.6 6.6
Columnea ciliata Jun 10.0 18.5 10.1 10.2 10.6 10.8 10.8 9.5
Columnea ciliata Mar 66.7 16.8 67.2 67.2 64.6 66.0 66.6 69.1
Columnea ciliata May 5.6 18.4 5.9 6.5 5.2 4.8 5.2 4.9
Columnea ciliata Nov 0.0 18.6 1.7 4.8 1.2 1.9 0.5 3.9
Columnea ciliata Oct 0.0 18.6 0.5 2.3 1.3 1.4 0.8 2.8
Columnea ciliata Sep 0.0 18.2 1.3 1.7 2.0 2.6 1.7 1.2
Columnea kucyniakii Apr 6.7 9.3 7.3 7.8 6.8 6.4 7.3 8.2
Columnea kucyniakii Aug 0.0 9.1 1.0 0.7 1.3 1.2 1.0 1.7
Columnea kucyniakii Dec 12.5 9.6 12.7 16.7 13.2 14.6 10.7 15.2
Columnea kucyniakii Feb 43.8 9.8 43.7 42.8 43.5 42.9 42.6 41.4
Columnea kucyniakii Jan 25.0 9.4 24.7 23.0 25.4 24.0 24.0 25.1
Columnea kucyniakii Jul 0.0 10.1 1.3 1.6 1.5 1.9 0.5 1.6
Columnea kucyniakii Jun 0.0 9.3 1.0 1.3 1.4 0.8 1.7 0.5
Columnea kucyniakii Mar 16.7 10.5 16.4 18.0 15.9 15.3 17.3 18.5
Columnea kucyniakii May 0.0 9.9 0.6 1.7 0.8 1.6 1.1 1.2
Columnea kucyniakii Nov 11.1 8.6 11.1 14.6 11.1 13.6 10.0 14.6
Columnea kucyniakii Oct 0.0 9.4 1.1 2.1 1.6 2.2 0.9 2.3
Columnea kucyniakii Sep 0.0 9.7 1.2 2.5 0.6 1.0 1.4 1.7
Columnea medicinalis Apr 6.7 18.2 6.1 5.4 7.1 7.6 7.2 9.9
Columnea medicinalis Aug 25.0 19.1 26.5 25.6 25.8 23.2 23.3 24.8
Columnea medicinalis Dec 0.0 18.6 0.8 3.6 0.7 4.4 1.6 4.7
Columnea medicinalis Feb 12.5 18.8 12.6 12.1 13.1 13.1 11.6 10.5
Columnea medicinalis Jan 16.7 19.5 16.8 15.4 17.1 17.2 16.2 18.2
Columnea medicinalis Jul 33.3 18.8 33.1 33.1 33.8 33.5 32.4 33.9
Columnea medicinalis Jun 30.0 18.2 29.3 29.2 30.3 29.9 30.5 30.6
Columnea medicinalis Mar 25.0 19.1 25.3 24.5 25.4 26.4 24.3 26.6
Columnea medicinalis May 11.1 18.8 10.8 10.6 11.4 11.6 11.2 10.6
Columnea medicinalis Nov 11.1 19.0 11.3 13.5 11.4 13.3 11.1 12.8
Columnea medicinalis Oct 28.6 18.9 29.4 28.9 27.5 29.2 28.0 31.0
Columnea medicinalis Sep 6.7 19.3 6.5 8.6 5.3 7.5 6.2 9.6
Columnea picta Apr 13.3 16.2 12.9 13.1 13.8 15.6 15.1 13.9
Columnea picta Aug 8.3 17.0 7.1 10.0 9.0 7.2 8.5 7.8
Columnea picta Dec 12.5 16.5 12.2 13.9 12.0 16.2 12.7 15.9
Columnea picta Feb 12.5 17.1 12.7 13.2 13.4 12.7 11.9 12.2
Columnea picta Jan 25.0 16.6 25.5 23.1 26.0 22.9 24.1 24.0
Columnea picta Jul 13.3 16.4 14.4 13.2 12.1 13.2 12.9 14.3
Columnea picta Jun 35.0 17.1 35.1 34.4 34.0 34.8 33.8 34.1
Columnea picta Mar 0.0 16.8 0.4 3.1 1.5 3.6 2.1 2.5
Columnea picta May 33.3 17.1 32.6 33.3 32.9 31.6 34.2 33.2
Columnea picta Nov 0.0 17.0 1.9 2.8 1.5 3.4 1.2 4.3
Columnea picta Oct 0.0 17.2 0.9 2.1 1.8 3.2 1.1 2.1
Columnea picta Sep 13.3 16.4 13.4 12.8 13.2 13.0 13.2 14.9
Columnea strigosa Apr 0.0 12.8 1.3 2.0 1.2 1.6 0.5 1.7
Columnea strigosa Aug 16.7 14.1 14.9 16.0 15.9 14.8 17.9 16.1
Columnea strigosa Dec 0.0 13.6 1.5 3.9 1.1 4.5 1.7 3.4
Columnea strigosa Feb 12.5 13.7 13.4 12.5 12.1 11.4 13.1 11.9
Columnea strigosa Jan 8.3 13.3 10.1 8.3 7.2 6.0 8.4 5.6
Columnea strigosa Jul 26.7 13.0 26.1 26.0 26.7 26.4 26.8 27.2
Columnea strigosa Jun 40.0 13.2 40.2 37.4 39.6 39.8 39.9 39.9
Columnea strigosa Mar 0.0 12.7 0.6 2.4 1.6 3.9 1.1 2.2
Columnea strigosa May 22.2 14.0 23.0 23.1 23.3 21.7 21.3 20.6
Columnea strigosa Nov 0.0 13.0 0.6 3.5 1.8 4.4 1.7 3.4
Columnea strigosa Oct 0.0 12.8 1.2 2.5 1.7 3.5 1.1 3.0
Columnea strigosa Sep 0.0 13.6 1.1 2.5 1.5 1.9 0.7 2.5
Drymonia collegarum Apr 13.3 13.8 14.9 11.8 13.7 13.8 13.4 13.4
Drymonia collegarum Aug 16.7 14.6 17.5 16.3 17.4 15.0 17.6 15.3
Drymonia collegarum Dec 25.0 15.1 25.2 26.4 25.2 29.8 24.6 28.3
Drymonia collegarum Feb 6.2 13.9 6.6 6.5 6.6 6.6 6.2 6.8
Drymonia collegarum Jan 16.7 14.5 17.4 15.1 17.5 16.4 16.2 16.6
Drymonia collegarum Jul 13.3 13.6 15.6 14.0 13.6 14.9 13.6 13.2
Drymonia collegarum Jun 15.0 14.2 16.0 14.2 13.8 15.1 16.6 14.2
Drymonia collegarum Mar 8.3 14.9 7.5 9.1 9.9 9.4 8.6 9.6
Drymonia collegarum May 22.2 14.4 22.0 22.4 21.5 21.8 21.6 22.8
Drymonia collegarum Nov 11.1 14.9 10.2 15.4 12.8 13.9 12.3 14.3
Drymonia collegarum Oct 7.1 14.2 7.6 8.9 7.8 7.9 6.3 8.8
Drymonia collegarum Sep 13.3 14.6 13.5 14.1 13.3 13.5 13.7 13.0
Drymonia tenuis Apr 13.3 15.1 14.2 13.4 12.7 13.9 14.0 12.3
Drymonia tenuis Aug 8.3 15.0 8.4 7.9 9.5 7.1 8.5 9.0
Drymonia tenuis Dec 0.0 14.8 2.1 6.2 1.5 4.2 3.1 4.5
Drymonia tenuis Feb 18.8 15.6 19.5 19.8 20.3 17.6 18.0 19.4
Drymonia tenuis Jan 16.7 15.1 16.4 15.8 16.6 15.1 18.9 14.8
Drymonia tenuis Jul 20.0 14.9 20.6 19.3 21.4 20.3 20.1 21.3
Drymonia tenuis Jun 20.0 15.2 19.7 20.3 20.8 19.4 19.8 20.5
Drymonia tenuis Mar 25.0 14.8 26.3 26.6 23.7 28.1 24.6 27.2
Drymonia tenuis May 22.2 15.3 22.2 22.9 23.2 22.6 22.3 23.1
Drymonia tenuis Nov 0.0 15.2 1.8 4.1 2.0 2.9 1.9 0.9
Drymonia tenuis Oct 7.1 14.7 8.2 9.6 7.4 7.2 6.4 7.2
Drymonia tenuis Sep 6.7 15.3 7.0 7.8 7.4 8.2 6.7 5.9
Drymonia teuscheri Apr 26.7 17.4 28.0 27.3 24.9 27.5 27.6 27.3
Drymonia teuscheri Aug 25.0 17.3 24.7 26.0 26.2 21.8 23.6 23.9
Drymonia teuscheri Dec 0.0 17.3 0.8 4.6 1.7 5.7 2.1 2.2
Drymonia teuscheri Feb 6.2 17.8 7.5 7.2 5.6 6.7 5.5 6.7
Drymonia teuscheri Jan 16.7 17.6 16.3 14.8 16.5 15.5 16.8 16.8
Drymonia teuscheri Jul 13.3 17.1 13.4 14.3 11.4 14.8 14.1 13.4
Drymonia teuscheri Jun 20.0 18.0 19.9 20.4 20.6 20.7 19.9 20.9
Drymonia teuscheri Mar 25.0 17.9 27.0 26.4 24.2 25.8 24.9 25.8
Drymonia teuscheri May 33.3 18.2 32.9 32.4 33.0 31.0 33.7 32.7
Drymonia teuscheri Nov 11.1 17.4 11.1 13.4 10.4 13.3 12.4 14.9
Drymonia teuscheri Oct 14.3 17.3 14.6 15.7 13.8 16.1 14.9 15.6
Drymonia teuscheri Sep 6.7 17.8 7.0 7.4 6.7 7.9 6.8 8.2
Gasteranthus lateralis Apr 13.3 14.6 12.4 12.1 12.8 13.7 13.2 12.3
Gasteranthus lateralis Aug 16.7 15.6 18.2 16.1 16.9 14.9 16.4 16.4
Gasteranthus lateralis Dec 25.0 15.3 26.2 27.1 25.2 28.6 23.6 27.4
Gasteranthus lateralis Feb 12.5 15.3 12.9 11.7 13.2 11.5 12.5 12.9
Gasteranthus lateralis Jan 33.3 15.7 34.0 32.8 34.0 31.7 33.8 33.3
Gasteranthus lateralis Jul 0.0 14.7 2.1 1.5 1.3 1.2 1.5 1.9
Gasteranthus lateralis Jun 5.0 15.3 5.1 5.5 5.2 4.2 6.6 5.3
Gasteranthus lateralis Mar 8.3 15.8 8.4 8.4 8.7 10.0 8.3 8.7
Gasteranthus lateralis May 5.6 15.2 6.5 6.1 5.7 4.9 7.7 5.7
Gasteranthus lateralis Nov 44.4 15.9 44.2 44.5 44.2 46.3 43.4 47.3
Gasteranthus lateralis Oct 28.6 15.7 28.1 27.7 27.3 29.8 27.8 29.0
Gasteranthus lateralis Sep 13.3 14.9 13.8 14.1 12.4 12.4 13.4 13.5
Gasteranthus quitensis Apr 6.7 13.5 7.8 6.2 6.3 8.9 6.8 6.3
Gasteranthus quitensis Aug 16.7 12.4 15.9 16.9 18.0 14.6 15.8 15.5
Gasteranthus quitensis Dec 12.5 12.2 11.9 17.5 12.3 17.3 13.4 16.0
Gasteranthus quitensis Feb 18.8 12.2 19.7 16.8 19.4 18.7 19.0 20.5
Gasteranthus quitensis Jan 41.7 13.1 42.4 40.1 42.1 39.9 42.0 39.7
Gasteranthus quitensis Jul 6.7 12.0 6.6 7.7 7.1 8.4 6.5 8.5
Gasteranthus quitensis Jun 0.0 12.5 1.2 1.5 0.7 0.6 1.2 1.8
Gasteranthus quitensis Mar 8.3 12.2 6.7 8.7 7.8 9.4 8.1 10.6
Gasteranthus quitensis May 0.0 13.2 1.2 1.2 0.5 1.0 1.5 1.3
Gasteranthus quitensis Nov 11.1 13.1 11.6 14.7 12.1 13.4 13.1 12.1
Gasteranthus quitensis Oct 21.4 11.8 22.1 23.3 21.0 22.7 20.6 21.7
Gasteranthus quitensis Sep 13.3 12.6 13.6 14.5 12.0 13.5 13.0 12.5
Glossoloma oblongicalyx Apr 0.0 15.4 1.4 1.4 0.9 1.2 0.7 1.1
Glossoloma oblongicalyx Aug 33.3 14.8 32.4 32.5 30.6 32.1 33.8 32.5
Glossoloma oblongicalyx Dec 0.0 15.9 1.7 3.4 1.8 3.8 0.8 4.1
Glossoloma oblongicalyx Feb 0.0 15.9 0.7 1.0 1.3 1.8 0.8 0.4
Glossoloma oblongicalyx Jan 0.0 15.8 0.7 1.2 1.9 0.8 1.4 0.6
Glossoloma oblongicalyx Jul 33.3 16.1 31.5 33.1 32.7 34.0 35.2 33.9
Glossoloma oblongicalyx Jun 10.0 16.0 10.4 8.8 9.7 10.8 10.2 9.1
Glossoloma oblongicalyx Mar 0.0 15.2 1.6 3.3 0.8 2.4 0.6 1.1
Glossoloma oblongicalyx May 11.1 15.3 11.3 11.3 11.5 9.8 10.5 11.1
Glossoloma oblongicalyx Nov 22.2 16.0 21.0 21.6 22.1 24.5 21.7 24.2
Glossoloma oblongicalyx Oct 35.7 15.7 37.2 36.3 33.6 37.5 36.1 35.3
Glossoloma oblongicalyx Sep 33.3 15.3 32.7 34.1 32.3 32.4 34.2 36.3
Glossoloma purpureum Apr 6.7 13.2 8.3 8.0 8.2 8.7 6.9 7.9
Glossoloma purpureum Aug 16.7 12.5 17.7 16.6 15.5 14.4 17.5 15.2
Glossoloma purpureum Dec 12.5 12.2 12.9 18.6 12.6 15.6 13.9 15.9
Glossoloma purpureum Feb 12.5 12.3 11.5 12.1 11.8 12.3 11.6 12.3
Glossoloma purpureum Jan 16.7 12.6 17.6 15.6 17.6 15.3 15.1 15.8
Glossoloma purpureum Jul 20.0 12.4 20.0 20.4 20.1 20.6 19.9 20.6
Glossoloma purpureum Jun 10.0 12.7 10.6 10.4 10.7 9.8 11.1 8.8
Glossoloma purpureum Mar 0.0 13.0 1.0 2.2 1.2 2.3 0.8 2.9
Glossoloma purpureum May 5.6 13.0 4.9 7.2 7.5 5.7 6.3 7.0
Glossoloma purpureum Nov 11.1 13.1 11.4 13.2 9.5 14.5 13.0 11.8
Glossoloma purpureum Oct 14.3 12.3 15.7 16.7 13.4 15.0 14.2 17.5
Glossoloma purpureum Sep 26.7 13.2 27.5 27.3 26.8 27.1 26.8 29.8
Kohleria affinis Apr 46.7 17.3 46.7 46.2 45.3 49.1 46.5 50.0
Kohleria affinis Aug 8.3 17.5 8.7 9.0 6.9 7.7 9.0 10.1
Kohleria affinis Dec 12.5 16.6 11.2 13.7 14.3 17.5 10.5 16.6
Kohleria affinis Feb 25.0 17.6 24.0 23.9 24.0 26.2 25.3 25.1
Kohleria affinis Jan 8.3 18.4 10.0 7.9 8.3 9.8 9.7 6.7
Kohleria affinis Jul 0.0 17.4 1.5 2.6 1.1 0.9 0.4 1.5
Kohleria affinis Jun 15.0 17.6 15.0 14.4 14.2 14.2 15.5 14.3
Kohleria affinis Mar 41.7 17.5 41.5 40.7 38.5 45.4 40.5 42.8
Kohleria affinis May 27.8 18.2 29.4 26.8 27.3 28.2 28.3 28.6
Kohleria affinis Nov 0.0 17.9 2.2 2.1 2.0 3.4 2.1 4.4
Kohleria affinis Oct 0.0 17.1 0.5 2.3 1.9 1.4 1.0 1.6
Kohleria affinis Sep 6.7 17.1 7.3 7.0 6.6 6.9 7.1 7.4

20 Flowering correlations in predictions